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ABSTRACT 

We present higk-resoiution X-ray spectra of cloud-shock interaction regions 
in the eastern and northern rims of the Galactic supernova remnant Puppis A, 
using the Reflection Grating Spectrometer onboard the XMM-Newton satellite. 
A number of emission lines including Ka triplets of He-like N, O, and Ne are 
clearly resolved for the first time. Intensity ratios of forbidden to resonance 
lines in the triplets are found to be higher than predictions by thermal emission 
models having plausible plasma parameters. The anomalous line ratios cannot 
be reproduced by effects of resonance scattering, recombination, or inner-shell 
ionization processes, but could be explained by charge-exchange emission that 
should arise at interfaces between the cold/warm clouds and the hot plasma. 
Our observations thus provide observational support for charge-exchange X-ray 
emission in supernova remnants. 

Subject headings: atomic processes — ISM: abundances ISM: individual ob- 
jects: Puppis A — ISM: supernova remnants — X-rays: ISM 
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Here, we present high-resolution X-ray spectra of the BEK/NK features using XMM- 
Newton!^ Reflection Grating Spectrometer (RGS: den Herder et al. 2001). From the RGS 
spectra, forbidden-to-resonance line ratios in Hea triplets are found to be anomalously en- 
hanced, especially at the BEK. We show that this anomaly could be naturally caused by the 
presence of CX emission. 


2. Observations and Spectral Analysis 

The NK and BEK were observed by XMM-Newton on 2003 April 17 (Obs.ID 0150150101) 
and 2003 May 21 (Obs.IDs 0150150201 and 0150150301), respectively, in order to obtain 
high-resolution spectra using the RGS instrument. The dispersion directions of the RGS are 
76°. 4 (NK) and 132°.3 (BEK) counterclockwise from the north, as shown in Fig. 1, which 
is an X-ray image of Puppis A generated from existing XMM-Newton and Chandra data. 
In this paper, we mainly focus on the RGS data, while we also utilize data taken by the 
European Photon Imaging Camera (EPIC: Turner et al. 2001; Striider et al. 2001) to support 
our RGS analyses. The exposure times after removing time periods affected by soft protons 
are 13.5 ks and 20.8 ks for the NK and BEK, respectively. All the raw data are processed 
using version 11.0.0 of the XMM Science Analysis Software and the latest calibration data 
files , available at the time of the analysis. 

Since the RGS is a slitless spectrometer, off-axis sources along the dispersion direction 
are detected at wavelength positions shifted with respect to the on-axis source. Spatial 
displacement of V corresponds to a wavelength shift of 0.138A(or 4eV/arcmin at 0.6 keV) 
for the first spectral order. Because the BEK and the NK are fairly compact (2'-3') features 
with much brighter surface brightness compared to their surroundings, the RGS is capable 
of producing high-resolution spectra for them, with an order-of-magnitude better resolution 
than nondispersive CCDs. In fact, there are a number of successful RGS observations not 
only of moderately extended SNRs in LMC/SMC (e.g., Rasmussen et al. 2001) but also of 
a bright knot in the northwestern rim of SN 1006 (Vink et al. 2003). 

As shown in Fig. 1, we divide the eastern RGS field of view into four sectors spaced 
by O'. 8 along the cross-dispersion axis, while we extract one spectrum from the l'.6-width 
region for the north. In this way, we mitigate the photon number difference between the 
BEK and the NK. We smooth RGS responses originally designed for point sources, based 
on emission profiles along the dispersion direction with the rgsrmf smooth software. For the 
input to the software, we arrange the X-ray images such that Chandra covers regions around 
the BEK/NK while XMM-Newton covers the remaining regions, and we take account of 
vignetting effects of XMM-Newton’s X-ray telescope. We generate energy-dependent RGS 
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the normalization ratios of 0 Ly7/0 Ly/3 and 0 Lyd/O Ly/3 to 0.3 and 0.1 as are expected 
at kT e ~0.5keV (Hwang et al. 2005). We note that the kT e values of the bremsstrahlung 
component are mainly determined by continuum emission above 1.5 keV. Thus, fixing ratios 
of Gaussian normalizations does not affect the kT e measurements. Line centers of the most 
prominent 21 Gaussians are left as free parameters, while those for C He/3 and two Fe L lines 
(G 1 and D) axe systematically shifted with respect to those of their neighboring lines (H and 
E). Other weak Gaussians’ centers are fixed to the theoretically expected values (Smith et 
al. 2001). The widths of all Gaussians are fixed to zero, since significant broadening is not 
required from a statistical point of view. With this fitting strategy, we obtain fairly good 
fits for all the spectra as shown in Fig. 2. The fit results are summarized in Table 1, where 
we omit weak Gaussians. 

The values of jVh are consistent with recent X-ray measurements (Hwang et al. 2005, 
2008; Katsuda et al. 2010). The kT c values are somewhat lower than previous results. 
However, prominent line intensities, which are essential for the discussion below, are not 
affected by the temperature difference, as the fraction of underlying continuum is very small. 
Line intensity ratios are generally reproduced by thermal emission models having plausible 
plasma parameters. On the other hand, the ratios of forbidden-to-resonance (F/R) lines in 
Hea triplets are apparently inconsistent with expectations of thermal emission models, as 
shown in Table 2. This is particularly evident for O ions in the BEK. To see if the forbidden 
line is stronger or the resonance fine is weaker with respect to the thermal prediction inferred 
by other lines and continuum emission, we fit the X-ray spectra (in the E2 region) with an 
absorbed vpshock model (Borkowski et al. 2001), excluding the energy band of O Her* 
triplets. After fitting, we recover the triplet lines and find that either the intensity of the 
forbidden line inferred from the vpshock model is weaker than the data or the intensity of 
the resonance line is stronger, or both of them. 


3. Discussion 

We have presented high-resolution X-ray spectra of the BEK/NK features in the Galactic 
SNR Puppis A. The forbidden and resonance fines in the Hea triplets are clearly resolved, 
and their intensity ratios (F /R) are found to be generally higher than thermal predictions. 
This anomaly is particularly evident in the BEK, while the RGS spectrum of the NK as well 
as Einstein’s FPCS spectra of the northeastern portion (Winkler et al. 1981a, b) are more 
close to thermal predictions. 


1 We use the same labeling for Fe XVII L lines as in Gillaspy et al. (2011). 
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Yamaguchi et al. 2009; Uchida et al. 2012, and references therein). The recombining plasmas 
axe characterized by strong radiative recombination continua (RRC) and enhanced Lya/Hea 
ratios. We find that signatures of RR.Cs (i.e., strong recombination edges) are not evident in 
the RGS and MOS spectra in Fig. 2, although the Fe L and/or other fine emission might make 
it difficult for us to detect such signatures. Also, the Lya/Hea ratios for O axe all around 0.6, 
which can be achieved by either under-ionized or equilibrium plasmas when kT e > 0.2keV 
(Smith et al. 2001). In addition, recombination processes would enhance the Fe L lines ratio 
of (F+G+H)/(C+D-(-E) as high as >25 (Liedahl et al. 1990), which is inconsistent with 
our RGS measurements of 2.0±0.1 at the E2 region. Furthermore, our spectral fitting of the 
combined RGS and MOS spectra with a recombining plasma model (the cie model in SPEX: 
Kaastra, et al. 1996) failed to reproduce the entire X-ray spectrum; the model requires too 
low electron temperature to explain emission above 1 keV. These investigations indicate that 
the plasma here is not recombining. 

Signatures of inner-shell ionization processes of Li-like ions can be found as Li-like 
satellite lines as well as enhanced forbidden-to-intercombination (F/I) ratio compared with 
collisionally excited emission (e.g., Porquet et al. 2010). In the RGS spectra, there is no indi- 
cation of satellite lines, however. Also, the measured O VII F/I ratios in Table 2 marginally 
consistent with thermal expectations. These facts led us to consider that inner-shell ioniza- 
tion processes are not working efficiently in this region. We note that, whereas line emission 
from Li-like O is seen in the far-ultraviolet spectrum of the BEK (Blair, et al. 1995), the 
abundance of Li-like O would be small in the X-ray emitting region. 

We then assess the feasibility of the CX scenario by calculating the expected CX flux, 
following Lallement (2004). The volume emissivity of CX is expressed as Pcx = 0cx Ur nj V). 
We let <7cx he the CX cross section between neutral H and ions of interest, tir the neutral 
H density, n; the ion density, and V t the relative H-ion velocity. We first consider the O 
VII forbidden line, so that <tcx is 3.3 xlO -15 cm -2 (Bodewits, et al. 2007). The value of 
Ur can be taken from the density of cold/warm clouds immersed in X-ray-emitting plasma, 
which is ~50cm -3 (Teske & Petre 1987). The value of n\ is the density of O VIII in the 
hot plasma, which is ~2x 10 -4 cm -3 for the BEK (Hwang et al. 2008; Arendt et al. 2010). 
We assume that V T = 500 km s -1 , roughly the shock velocity that can produce the X-ray- 
emitting plasma in the BEK. These parameters give P C x ~ 1- 7 x 10 9 photons cm -3 s -1 . The 
CX-emitting volume is calculated by the thickness of the CX layer, lex, times the interface 
area between the clouds and the hot plasma. The value of lex is equated to the mean free 
path for H-proton CX. Thus, lex is of the order of l/cm p , with a being the H-proton CX 
cross section (10 -15 cm -2 : McClure 1966) and n p being the proton density in the plasma 
(4cm -3 : Arendt et al. 2010). Considering that only ~30% neutral H can charge-exchange 
before getting collisionally ionized (Lallement 2004) for the BEK’s plasma condition, we 
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M82. In this context, we conclude that CX emission is a promising mechanism for explaining 
the anomalous F/R ratio observed in Puppis A. 

One may expect to find other spectral CX signatures in addition to the Hea line ratios. 
However, CX spectral properties axe strongly dependent on V T and target neutrals (e.g., 
Beiersdorfer, et al. 2001, 2003), and investigation of CX emission is still ongoing. Therefore, 
further discussions call for more sophisticated CX emission modeling, which is beyond the 
scope of this paper and is left as future work. We also expect future observations with 
the non-dispersive Soft X-ray Spectrometer (SXS: Mitsuda et al. 2010) onboard the Astro- 
H satellite (Takahashi et al. 2010) to provide additional interesting information especially 
regarding spatial structures. 


4. Conclusion 

High-resolution X-ray spectra of the cloud-shock interaction regions, the BEK and the 
NK, in Puppis A revealed anomalous Hea triplet ratios: the O Hea F/R line ratios are found 
to be ~2 in the BEK. This anomalous ratio could be naturally interpreted as a result of the 
presence of CX emission, while resonance-scattering effects might be non-negligible as well. 
CX X-ray emission has been recently found in or proposed for many astrophysical sites such 
as comets (e.g., Lisse et al. 1996), diffuse soft X-ray background (e.g., Fujimoto et al. 2007; 
Ezoe et al. 2011), star-forming regions (Townsley et al. 2011a, b), and starburst galaxies (e.g., 
Tsuru et al. 2007; Ranalli et al. 2008; Konami et al. 2011; Liu et al. 2012). The Astro-Ifs 
SXS will allow us to obtain high-resolution, spatially-resolved spectra of many more objects, 
shedding additional light on the role of CX processes in X-ray astrophysics. 
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experiments of CX X-ray emission. S.K. and S.Y. are supported by the Special Postdoc- 
toral Researchers Program in RIKEN. This work is partly supported by a Grant-in- Aid for 
Scientific Research by the Ministry of Education, Culture, Sports, Science and Technology 
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Table 1. Spectral-fit parameters 


Component 

Parameter 



Region 





El 

E2 

E3 

E4 

N 

Absorption 

ATh (10 i4l cm“ a ) 

2.70+0.01 

280+0.01 

2.80+0.01 

2.85±0.01 

2.58+0.01 

Bremsstrahlung 

kT 9 (keV) 

0.32 ±0.01 

0.3i±0.01 

0.33*0.01 

0.32*0.01 

0.30*0.01 


Normalization® 

5264±S3 

9124±133 

6744±117 

7808±117 

5397+90 

Gaussian: C Lya+Sl L-f-S L 

Center (eV) 

390.8t£l 


366.71®;! 

365.412 * 

387.412! 


Normalization 0 

822+80 

I486* 118 

1257±139 

1723±139 

380+52 

N Hea (f) 

Center (eV) 


419.2±0.3 

418.4+0.5 

419.7+0.5 

420.6+”;! 


Normalisation 0 

100+18 

201 ±26 

145±28 

243±33 

82+16 

N Hea (r) 

Center (eV) 

429.71“-® 

43o.i12;| 

429.9±0.5 

429.7±0.5 

430.6+0.5 


Normalization 0 

146+17 

223±24 

1 69 ±23 

199±27 

72+15 

c Lyfl 

Center (eV) 

•*34.91“-’; 

436.1 12" 5 

435.212;! 

*34.71“;! 

*35.21“;! 


Normalization® 

60±14 

87±18 

65±18 

87±22 

51+13 

C Ly-r 

Center (eV) 

459® 

459® 

459® 

459® 

459® 


Normalization® 

25+8 

31±10 

24±10 

14+12 

< 8 

N Lya-hHe/J 

Center (eV) 

499.4+0.2 

500.212.J 

499.812! 

499.6l“J 

499.612^ 


Normalization® 

143±9 

254±13 

222±12 

252+14 

162+11 

O Hea (f) 

Center (eV) 

559.8+“* 

560.7±0.1 

560.5±0.G 

560.2±0.5 

561.2+0.2 


Normalization® 

705+27 

1167±36 

1042±36 

1119 + 36 

417+23 

O Hea (1) 

Center (eV) 

686.81^ 

D68.3l5;2 

570.71“;! 

668.512* 

567.51“^ 


Normalization® 

107±24 

190±30 

298±31 

265±34 

168+22 

0 Hea (r) 

Center (eV) 

573.1±0.2 

574.i12;| 

573.4±0.1 

572.9±0.1 

574.3+2;® 


Normalization® 

601 ±25 

771±30 

527±31 

646+33 

492+21 

O Lya 

Center (eV) 

662.7l^J 

653.4+0.2 

603.0+0.2 

652.8+0.2 

653.4+0.4 


Normalization® 

97»±38 

1223±47 

1154 + 4 0 

1329+51 

631+36 

O He/? 

Center (eV) 

664.4+0.8 

685.012J 

665.4+2! 

666.0+}'! 

667.0+1 


Normalization® 

161±24 

285 ±29 

197+27 

136+29 

110+17 

O He-r 

Center (eV) 

699.01JJ 

700.7±0.8 

698.91** 

698.21!;! 

700.7+1.1 


Normalization® 

68+8 

74±10 

65±U 

74+13 

42+9 

Fe L (G+H) 

Center c (eV) 

724.9l“;J 

725.1+“;! 

724.9+2;! 

725.2+1* 

725.2+“;! 


Normalization® 


3621®* 

33lliS 

332 IS 

issliT 

Fe L (F) 

Center (eV) 

739.312! 

738.51“* 

739.3l“! 

738.11“;! 

738.912;! 


Normalization® 

89±11 

156±13 

94±12 

99+12 

65+8 

O Ly/5 

Center (eV) 

773.612;! 

773.8+0.3 

773.312;! 

772.912! 

773.7+0.5 


Normalization® 

151±7 

187±8 

160±7 

169+8 

97+6 

Fe L (D+E) 

Center^ (eV) 

810.6±0.8 

810.1±0.2 

810.0±0.4 

810.0+0.4 

810.1+0.7 

Normalization® 

741!° 

noiii 

98+ 11 

QQ't H 
yu ~io 


Fe L (C) 

Center (eV) 

828.61“-* 

826.1±0.5 

825.6±0.5 

826.1+0.5 

826.2+“;® 


Normalization® 

105+6 

X73±0.8 

135±7 

152+8 

78+5 

Ne Hea (f)+Fe L 

Center (eV) 

902. 0± 0.6 

903.312^ 

902.612J 

902 612;! 

903.5+S ? 


Normalization® 

122±6 

185±7 

156±6 

175+7 

75+4 

Ne Hea (r)+Fe L 

Center (eV) 

921.212;! 

921.612! 

921.21“^ 

921.8+2"! 

922.4+0.2 


Normalization® 

141±0 

224±7 

168±8 

188±7 

111+5 

Ne Lya+Pe L 

Center (eV) 

1023.6+*"! 

1022.8+!'! 

1020.5+2"! 

1022.1+“! 

1020.41!;* 


N orm alization ° 

83±6 

97±6 

96 ±5 

102+6 

35+3 

x*/d.o.f. 


1423/968 

1740/1151 

1698/1114 

1684/1181 

949/612 


Noto. — ®In units of 10 4 photons cm -2 a” 1 . 

^ Fixed values. 

c Line centers of K and E are shown, while those of G and D are systematically shifted by +2eV from H and E, respectively. 
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Fig. 2. — XMM-Newton spectra of the five regions in Fig. 1. The data covering below 0.6 keV, 
0.65-1.5keV, and above 1.1 keV are the first order RGS1+2, the second order RGS1+2, and 
MOS1+2, respectively. These data are simultaneously fitted with a phenomenological model 
(see text), and the best-fit models are shown in red. Lower panels show the residuals. The 
insets show close-up spectra for 0 Hea triplets with individual best-fit model components. 


